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Much progress has been made in the field of quantum computing using continuous variables 
over the last couple of years. This includes the generation of extremely large entangled cluster 
states (10,000 modes, in fact) as well as a fault tolerant architecture. This has led to the point 
that continuous-variable quantum computing can indeed be thought of as a viable alternative for 
universal quantum computing. With that in mind, we present a new algorithm for continuous- 
variable quantum computers which gives an exponential speedup over the best known classical 
methods. Specifically, this relates to efficiently calculating the scattering amplitudes in scalar bosonic 
quantum field theory, a problem that is believed to be hard using a classical computer. Building on 
this, we give an experimental implementation based on cluster states that is feasible with today’s 
technology. 

PACS numbers: 42.50.Ex, 03.70.+k, 42.50.Dv, 03.67.Lx 


I. INTRODUCTION 

For more than a decade now, continuous-variable (CV) 
quantum information mm has been a prominent sub¬ 
strate in implementing quantum technologies. Primar¬ 
ily this can be attributed to its largely Gaussian nature 
which invites simple and convenient mathematical cal¬ 
culations, as well as accessible experimental demonstra¬ 
tions, which are often deterministic in nature. Further¬ 
more, one can also use CVs as a key element in another 
promising architecture, known as hybrid quantum infor¬ 
mation [3j. 

The field of quantum computing |3J using CVs HO has 
also progressed significantly in the last few years. From 
its original conception in 1999 [5j, progress began to ac¬ 
celerate after a cluster state [6! version was established 
in 2006 [7, B], leading to something significantly more 
tangible for experimentalists. This resulted in numer¬ 
ous proof-of-principle demonstrations [SHE], currently 
culminating in a 10,000 node cluster m created ‘on- 
the-go’ along with a 60 node cluster created simultane¬ 
ously El- From a theoretical perspective, much progress 
has been made ]15U23j . including recently, an important 
fault tolerant architecture [ 23 , achieved by leveraging the 
Gottesman-Kitaev-Preskill (GKP) encoding [23]. How¬ 
ever, one area that is significantly underdeveloped is that 
of algorithms for a CV quantum computer. Thus far 
there only exists CV versions of quantum searching [26] 
and the Deutsch-Jozsa algorithm [271 - 130] . 

In this paper, we present an algorithm that simu¬ 
lates m the scattering amplitudes in scalar bosonic 
quantum field theory (QFT) using a continuous-variable 
quantum computer. In fact, we show one can obtain an 
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exponential speedup over the best known classical algo¬ 
rithms. A discrete version of the algorithm was originally 
shown in Refs. [321133] for a quantum computer based on 
qubits. Further work extended this result to fermionic 
QFTs 23, as well as using wavelets for multi-scale sim¬ 
ulations [351 . 

Typically, q and p are the CVs spreading across all 
real numbers. To encode them in qubits, one needs a 
whole register of qubits at each point in space. However, 
with CVs, there is a 1-to-l mapping to qumodes (the 
CV equivalent of a qubit). In fact it is arguable that a 
CV quantum computer is the natural choice for such a 
QFT problem given that the fields are continuous vari¬ 
ables. Thus, the value of the field at a given point in 
space can be mapped onto a qumode naturally. If qubits 
are used, instead, the qumode needs to be replaced by a 
register of M qubits which only allows the field to take 
on 2 ai discrete values. Brennen et al. describe both pos¬ 
sibilities in Ref. [35], although they do not explain how 
to implement the quartic phase gate with CVs, which 
we do here. Furthermore, the quartic vertex in wavelets 
becomes very complicated. Implementing it would re¬ 
quire gates acting on more than two modes (resulting in 
logarithmic overhead in complexity). Another benefit to 
our approach is in the development of the initial cluster 
state. Here we show how to create the initial CV cluster 
state as well as suggesting an experimental implementa¬ 
tion based on standard linear optics. Furthermore, we 
also note that in the preparation of the initial state we 
see a slight improvement over the original qubit approach 
of Ref. [32]. There they require 0(N 2 - 376 ) gates to engi¬ 
neer the ground cluster state; whereas in our scheme, we 
require slightly less than that, specifically, 0(N 2 ) gates. 

Our paper is structured in the following way. In Sec.[TTl 
we discretize space for a one-dimensional scalar bosonic 
QFT while leaving the field and time as continuous pa¬ 
rameters. Next, we show how to generate the initial 
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cluster state using only Gaussian operations in Sec. Ill 


In Sec. IV we outline the steps necessary to compute 
a scattering amplitude including the required measure¬ 
ment. We provide an explicit experimental implementa¬ 
tion in Sec. |Vj Finally, the benefits of our approach over 
classical methods are discussed in Sec. IVII 


II. DISCRETIZATION IN ONE-DIMENSION 


We consider a relativistic scalar field (j> in one spa¬ 
tial dimension including a quartic self-interaction. We 
shall outline the discretization specifically in the one¬ 
dimensional case so as not to clutter the notation un¬ 
necessarily, but generalization to higher dimensions is 
straightforward and is discussed in the supplementary 
material. We note that the field (f> is a function of x and 
t (time), <fi(x,t). All three parameters are continuous. In 
our approach, we discretize x , but not <fi or t. In the case 
of qubits, one would discretize x and 0, but not t. In 
classical lattice calculations, one discretizes all three </>, 
x , and t. 

In the continuum, the one-dimensional free scalar QFT 
is given by the Hamiltonian 




•2 ±2 
m cp 


(1) 


where <f> is the scalar field and 7r the conjugate momentum 
field. They obey commutation relations [^>(2;), 7r(a/)] = 
i8(x — x') where we choose units in which h = 1 . 

We discretize space by letting x = na , n = 
0,1,..., iV — 1, where a is the lattice spacing and L = Na 
is the finite length of the spatial dimension (L > o). We 
choose units in which a = 1, for simplicity, and denote 
Qn = 4>(x), P n = 7 r(x). The discretized variables obey 
standard commutator relations, [Q„,P m ] = i8 nm . The 
Hamiltonian becomes 


N-l 


Hn = 


P* + m 2 Q 2 n , 1 


N-l 


n—0 


+ o yi (Qn ~~ Qn+ if (2) 


n—0 


where we employed periodic boundary conditions and de¬ 
fined Q n = Q 0 . 

It is useful to define creation and annihilation opera¬ 
tors, A\ and A n , respectively, by A n = ( Q n + iP n )/y/2. 
They obey the commutation relations = S nrn 

and the Hamiltonian can then be written as 

H 0 = \ P t P+^Q t VQ (3) 

where P = [P 0 , Pi,.. ., Pn-i] t and Q = 

[Qo, Qi, ■ ■ ■ ,Qn-i] T ■ The eigenvalues of the matrix V 
and the components of the corresponding normalized 
eigenvectors e" are, respectively, w 2 = m 2 + 4 sin 2 
and ej! = _L e 27r *W iv ) = 0, 1. Notice that 

the massless case is special because it contains a zero 


mode (for m = 0, wo = 0), so the matrix V is not 
invertible. To avoid the problems that arise, we can shift 
the mass by a small amount ~ 1/IV, which vanishes in 
the continuum limit (N -A oo). 

We also wish to add a quartic interaction, Hi nt = 
dxcj ) 4 -A Q\ which necessitates the addition 

of a mass counter term H c , t . = Jo dxcj) 2 -A Yh n Qn 
due to renormalization, as explained in the supplemen¬ 
tary material. We find that for weak coupling, the phys¬ 
ically interesting case is stable for A > 0. 

To diagonalize the Hamiltonian, we introduce new cre¬ 
ation and annihilation operators, aj. and at, respectively, 
defined by a k = V^f( et Q)fc + 75Sy( etp )fc where e is 
the matrix of the eigenvectors. Notice that e is unitary, 
e^e = I. These operators obey standard commutation 
relations, [ak,aj] = Su and the free Hamiltonian reads 



In this form, it is straightforward to construct the states 
in the Hilbert space. 


III. INITIAL CLUSTER STATE PREPARATION 

For the initial cluster state, in Refs. [321 S3 the excited 
state was created after creating the ground state. This is 
difficult because it involves manipulating a large number 
of qubits. In our approach, we create a single photon 
state in a single mode before creating the cluster state. 
This is more accessible, as it involves creating the state 
|1) for a single mode. It can be done in a variety of 
ways, via a heralded single photon source, for instance. 
At the end of the computation, the held modes are all 
measured and the distribution of single photons across 
them determines the result. 

To begin with, we build the system with N oscillators 
representing the variables (Q n , Pn), n = 0,1,2,.... The 
nth oscillator has a Hilbert space constructed by suc¬ 
cessive application of the creation operator A' n on the 
vacuum |0) n , which is annihilated by A n . Here |0) n is 
shorthand for a product state of vacuum fields 

|0) = |0) o ®|0)i®---®|0)jv_i , (5) 

with A„|0) =0. For a scattering process, we are given an 
initial state typically consisting of a fixed number of par¬ 
ticles, usually two, which undergoes evolution and then 
a measurement is performed (detection of particles) on 
the final state. Both initial and final states asymptote to 
eigenstates of the free Hamiltonian Hq. Thus quantum 
computation starts with preparation of an eigenstate of 
H 0 . 

First, we consider the ground state of H 0 . It is the 
cluster state |fi) annihilated by all a k , he., afc|fl) =0 for 
k = 0,1,..., N — 1. It can be constructed from the vac¬ 
uum state ([5| by acting with the Gaussian unitary U\ 
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where a n = U'A n U. Noticing the relationship between 
the operators au and Ak we can use the Bloch-Messiah 
reduction [36] to determine U = VSW 1 as a decomposi¬ 
tion involving a multiport interferometer (V) followed by 
single mode squeezing ( S ) followed by a final multiport 
interferometer (W). These unitary operators can be re¬ 
alized with 0(N 2 ) gates |37| . This is in contrast to the 
qubit version [32] where they require 0(N 2 - 376 ) gates. 

To implement U we first perform the rotation 


iV-l 


Aq —> A' 0 = ^ A 

k =o 
N -1 

A' n = ^ 


k= 0 
N—l 


2irnk 

cos —4 


. . / v—v . 2'KTik 

An-u —► A N _ n = y ^ sin ———Ak 

k -o 


N 


( 6 ) 


where 1 < n < N/2, which can be expressed in terms 
of rotations each involving only a couple of oscillators. 
Notice that if N is even, Ajv /2 does not have a partner; we 
obtain A N / 2 —> ]G fc (— ) k Ak- Next, we squeeze each mode 

as A' n —> A" = cosh r n A' n + sinhr„A„ where e 2r " = oj n 
for n < N/2, and e _2r '* = ui n , for n > N/2. Finally, 
we untangle the pairs by rotating them, A 1 / —> ak where 
a 0 = A'/, a n = (. A" + iA%_ n )/V 2, and a N - n = ( iA " + 
A l /f_ n )/V2. Excited states can be constructed with the 
same number of gates, e.g., the single-particle state | k) = 
aj. | Q) can be constructed by acting upon the vacuum 
with A.j.. This turns the initial state of the fcth mode into 
a one-photon state, Aj,|0)fc, which can be accomplished 
in a variety of ways; see supplementary material. Having 
engineered Aj, |0)fe, we then apply the Gaussian unitary 
U\ to obtain the one-particle state 


a fc + |0) = E/t4| 0 ) 


(7) 


Extending the above to the engineering of multi-particle 
states, |fei, k 2 ,...) oc akjak 2 ' • ■ ■ 10), is straightforward. 


IV. QUANTUM COMPUTATION 

We wish to calculate a general scattering amplitude, 
which can be written as 

= (owt|Texp |*y dt(H int (t) + F c . t .(t))| \in) (8) 

in the limit T —> oo, where time evolution is defined with 
respect to the non-interacting Hamiltonian. 

We start by preparing the initial state | in) as in the 
previous section and define initial time as t = —T. Then 
we act successively with evolution operators of the form 


Time dependence is obtained via the free Hamiltonian, 
Qi(t) = e itHo Q t ( 0)e~ ltH ° (10) 

Therefore, the evolution can be implemented as 


U(t) = e itH ° e iSt ( H int+H 0 , t ,) e ~itH 0 


( 11 ) 


We deduce 

A = ( out\ 


e i6tH 0 e i6t(H irlt +H 0 , t .) 


N 


I in) (12) 


where we divided the time interval into N = ^ seg¬ 
ments. 

The coupling constants in ([9]) are turned on and off 
adiabatically. This is achieved by splitting the time in¬ 
terval \—T,T\ into three segments, [—T, — Ti], [—Ti,Ti], 
and [Ti, T], For t £ [—T, — Ti], we turn the coupling con¬ 
stants on by replacing A -A A (t), 6m -A 8m{t), so that 
A(— T) = Sm(—T ) = 0, and A(—Ti) = A, 6m{—Ti) = 6m. 
Then for t £ [—Ti,Ti] the coupling constants are held 
fixed. Finally, for t £ pi, T], they are turned off adiabat¬ 
ically by reversing the process in the first time interval. 
In the case of small A, the time dependence of the cou¬ 
pling constants can be chosen efficiently by making use 
of perturbative renormalization. Eqs. (A13) and (A14) 


inform the choice A(t) = A, 8m(t) log J4, for 


T-Ti' 

—T <t< — T\. 

The unitary operators e lStH ° and e lStHc t ‘ are Gaus¬ 
sian and can be implemented with second order nonlinear 
optical interactions and linear optics beam splitter net¬ 
works. The interaction is implemented through a quartic 
phase gate for each mode, 


AStHon 


n- 


,ilQn 


7 = 5t 


(13) 


The quartic phase gate may be implemented in a similar 
manner to the cubic phase gate previously proposed US- 
After evolution, we must project onto the state |ou£). 
This is similar to the state | in), and its construction de¬ 
pends on the number of desired particles. The latter are 
excitations created with aJ n , so in general, 


\°ut) = a/ ni a^ 2 


|n) = C7t«...|0) 


(14) 


It follows that the next step is to uncompute by applying 
the Gaussian unitary U (which is the inverse operation 
to the preparation of the initial state), and then measure 
the number of photons in each mode. The final uncom¬ 
pute step projects the set of output modes onto the Fock 
basis. Thus, the scattering amplitude calculation is a 
mapping from one set of field modes on the input to a 
separate set of field modes on the output, as expected. 
That is, for each click on the photodetector for mode n, 
there is an operator present in the final state (141. If 


U(t) = exp {i5t(H int (t) + H c , t ,{t))} 


(9) 


the QFT calculation involved an initial input state with 
two excitations spread across 100 field modes, say, then 
the entire calculation would involve two photons, for in¬ 
stance. We note that the calculation has made use of a 
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propagation phase shift/ squeezing 

interferometer 


FIG. 1: (Color Online) Sketch of an experimental setup for electromagnetic field modes used as qudits in a QFT 
calculation involving four field modes. The modes are encoded into electric field modes (colored red, blue, yellow, 
green), which are then prepared via beam splitters, swap gates, and squeezers for the compute stage. The compute 
stage consists of an interferometer, a quartic phase gate (black box, see Ref. [15]), and free propagation. An 
uncompute stage, which is the inverse of the preparation stage, and a detection stage in the Fock basis, yield the 
scattering amplitudes into the four QFT field modes. 


quartic phase gate up to this point, and thus technically 
speaking a non-Gaussian operation would not be neces¬ 
sary during this measurement step in order to achieve an 
exponential speedup over the classical QFT algorithm. 
However, in order to achieve high accuracy in the final 
result, photon number resolving detectors with high effi¬ 
ciency [38] would be desirable for the measurment phase. 

V. EXPERIMENTAL IMPLEMENTATION 

An example of the experimental implementation is 
given in Fig. |Tj For brevity the setup for calculating 
four space time points is given. For the electromagnetic 
field, the initial unitary rotation involves weighted beam 
splitters with the appropriate splitting to achieve the de¬ 
sired sums over the field operators (see appendix A.2). 
A swap gate is involved in the input state preparation 
stage. We note that a swap gate contains essentially the 
CV version of the CNOT operator along with parity op¬ 
erators [35], but in some cases the gate can be simplified 
to a beam splitter interaction PQ such as for the electro¬ 
magnetic field. Here we use a mode label swap operator, 
which is possible in systems with movable qubits, such 
as CV optical fields. Next, H c t . is quadratic in position 
quadrature operators, which can be implemented with a 
series of phase shifts pQ. The non-Gaussian piece of the 
computation is then the quartic phase gate contained in 
Hi n t, which can be implemented via repeated application 
of the photon number-dependent phase gate [15] . Lastly, 
the free propagation Hq can be implemented by a cali¬ 
brated free propagation before the uncompute stage. We 
note that the QFT field modes are encoded into the qu¬ 


dits which are themselves electromagnetic field modes, 
meaning that the free propagation contained in Hq is 
not arbitrary. It must conform to the calculated QFT 
free propagation distance, and phase stability must be 
maintained throughout. 


VI. CONCLUSION 

In conclusion, we developed a new algorithm for a 
continuous-variable quantum computer which gave an 
exponential speedup over the best known classical algo¬ 
rithms. This algorithm was the calculation of the scat¬ 
tering amplitudes in scalar bosonic quantum field theory, 
and as previously mentioned, arguably a natural choice 
for a continuous variable quantum computer to solve. At 
weak coupling, analytic calculations are possible, how¬ 
ever, at strong coupling no such calculations are generally 
available, and one has to rely on numerical techniques. 
A widely used framework is lattice field theory which is 
based on the discretization of space into a finite set of 
points. The complexity of classical computations on a 
lattice increases exponentially with the number of lattice 
sites pfO] . 

Quantum computations offer a distinct advantage (first 
shown in Ref. [33] for qubits, and here for qumodes), 
since complexity only grows polynomially. Using contin¬ 
uous variables we also see an advantage over the original 
qubit proposal; specifically, in the preparation of the ini¬ 
tial cluster state. There they required 0(N 2 - 376 ) gates 
to synthesize the ground state. However, in our scheme, 
we required slightly less, 0(N 2 ) gates. Finally, we also 
gave an example of an experimental implementation on a 
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continuous-variable cluster state quantum computer that 
calculated four space time points. We noted that such a 
scheme is feasible with current linear optical technology 
and consisted of a set of Gaussian operations along with 
the non-Gaussian quartic phase gate. 


The shift can be corrected by the addition of the counter 
term 

H c ,= 5 -f /W-►%£<& (AlO) 
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Appendix A: Renormalization 

Define the Green function G(t\,t 2 ) as 

G y (ti,t 2 ) = <0| T(Q i (ii)Q i (i2))|0>, (Al) 

where T denotes the time-ordering operator. It obeys 
K + V] G(h,t 2 ) = -H6(t i - t 2 ). (A2) 

Using the Fourier transform, 

G(t u t 2 ) = J (A3) 

we obtain 


with 5 m = — £ + G( A 2 ), i.e., the mass parameter in the 
Hamiltonian is not physical, but bare, 

ml = m 2 + 6 m = to 2 - -E £ — + 0{\ 2 ). (All) 
For large N, the sum can be approximated by an integral, 


£ = 



(A12) 


which has a logarithmic divergence at small m 2 (i.e., 
length scale 1/m large in units of lattice spacing, which 
is the physically interesting limit). We easily obtain 


£ = ^log 6 ^+0(m 2 ) 
o7T m z 


(A13) 


The bare mass is 


TOq = to 2 — £ + 0( A 2 ) = to 2 —— log Et + 0(A 2 , to 2 ) 

87t to 2 

(A14) 

Notice that for weak coupling (small A), the physically 
interesting case has TOg < 0 (a stable system, as long as 
A > 0). 


G(w) = i [-w 2 I + V] _1 = £ 


2 e « e n ) 


(A4) 


1. Ground State Construction 


exhibiting poles at w 2 = w 2 . 

When we switch on the interaction term, 


Hint a\ [ dx4> 4 ,| £ Q 


(A5) 


we have that at G(A) the Green function is corrected by 


SGij^h) = (0|T 


Qi(ti)Qj(t 2 ) / dtH int (t) 


For the Fourier transform, we obtain 

dco' 


| 0 >. 

(A6) 


<5G(w) = A[G(w)] 2 / UUTrG(w') (A7) 
./ 27T 


which leads to a shift of the poles, 
G(w) + <5G(w) = £ 

n 

where 


w 2 — w 2 - £ 


+ (A8) 


A f duj f \—a —i A v —a 1 . . . 

= 2N J 2^ E _ 2 = In E ^ ( } 


To find the required transformation U , we work as fol¬ 
lows. Notice that for n = 0, 


1 


cio — 


2 y/N 


N—l 

E 

k—0 


m H E) A fc + (yfm - ]= ) A\ 


(A15) 

where we used ujq = to. For n yf 0, we consider pairs 
( a n: a,N - n ). We have 


1 J 2nnk 
+ CL /V—r? = — 7= / cos ——— 

2 VN ^ N 


u n 4— i= ) Afc 


1 




Al 


i 27mA; 

On - a N - n = > sin 

fc=o 


2\/N 


N 


- VWn - 


At 


(A16) 


where we used w n = lun-u- 

The above expressions suggest that we transform A n 


into a n in three steps as detailed in Sec. Ill 
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2. Example: N = 4 

To illustrate the above algorithm, we consider the case 
in which space has been discretized to four points. The 
rotation (A' = OA) is described by the orthogonal ma¬ 
trix 



111 1 

y/2 0 -y/2 0 

1-11 -1 

0 V2 0 -y/2 


We have 


O — i ? 0 2 S 01 R 13 R-02 



(A17) 


(A18) 


where Rij(0) is a rotation in the ij-plane of angle 6 and 
Sij is the swap i j. Therefore the rotation O can be 
implemented with four two-mode unitaries. 

Next, we squeeze each mode as A' n —► A" = 
coshr n A' n + sinh r„A„ , where e 2r ° = ujq, e 2ri = uq, 
e 2r2 = W 2 , and e~ 2r3 = 0 J 3 . Notice that r 3 = —n, be¬ 
cause W 3 = U!-[. 

Finally, we perform the rotation, A" —► -^=(A" -MA 3 ), 
A 3 —> ~^(iA” + A 3 ), to arrive at the desired modes, 


n 

— 

n 

«2 = ^E( _1 ) 


“ 3 = 


cosh r 0 A n + sinh r 0 E A '~ 

n 

cosh r\A n + sinh rq E ■^n 

n 

coshr2A n + sinhr2 E ■^11 

n 

cosh 7*3 A n + sinhr 3 E 


(A19) 


Each of the above steps is implemented with a Gaus¬ 
sian unitary involving at most two modes. 


Appendix B: Excited States 

To generate the required one-photon state, two meth¬ 
ods can be used. One can first squeeze the vacuum of the 
kth mode with an optical parametric amplifier to 

S k (s) |0) fc , S k (a) = eH< 2 - A to (Bl) 

Then pass the squeezed state through a (highly trans¬ 
mitting) beam splitter of transmittance T, and place a 
photodetector on the auxiliary output port. A click of the 
detector heralds a successful photon subtraction, which 
is described by the non-unitary operator 


(B2) 


The transmittance has to be high so that the probability 
of detecting two or more photons is negligible. If no 
photon is detected, the process is repeated until a photon 
is detected. Finally, apply anti-squeezing Sl(s'). 

We obtain the state (unnormalized) 

St( s /) T AtA fe / 2 A fc 5 fc ( s )|o) fc (B3) 


If the squeezing parameters are chosen so that 


T = 


tanh s' 
tanh s 


(B4) 


then it is straightforward to show that (B31 is the desired 
state, 


Sl(s')T A i A ^ 2 A k S k (s)\0) k (x A{ |0) fc . (B5) 


Optionally, one may also use a heralded single pho¬ 
ton source. Such a source would consist of a parametric 
downconverter with a high efficiency heralding detector. 
To obtain exactly one photon when operating the source 
with high brightness (but on average less than one pair 
per pulse), the heralding detector would consist of a high 
efficiency photon number resolving detector, such as a 
transition edge sensor. 


Appendix C: Generalization to Arbitrary 
Dimensions 


Generalization to arbitrary spatial dimension d is 
straightforward. The free-scalar Hamiltonian in the con¬ 
tinuum reads 

H 0 = ^ J d d x [tt 2 + (V (/>) 2 + m 2 (j) 2 ] (Cl) 

where x G [ 0 ,L] d , with the fields obeying standard com¬ 
mutation relations, 


(x) , tt(x')] = iS d {x-x) 


(C2) 


Each coordinate Xi (i = 1,..., d) is discretized as before, 
Xi = Hid, Hi = 0,1,..., N — 1, and we define Q n = ^(x), 
P n = tt(x), A n = ^(Qn + iPn); where n G Z%. 

The Hamiltonian (Cl) can then be rendered in the 
form ([ 3 ]), where V has eigenvalues and corresponding 
normalized eigenvectors, 


ujI = to 2 + 4^]sin 2 -A 


e k N d/2 e 


(C3) 


where k G (the dual lattice). The eigenvectors 

form a unitary matrix. 

The discretized Hamiltonian is diagonalized as 


7?o = E Wk ( a k flk 


ker 


1 


y/1 -TT A * Ak ' 2 A k 


(C4) 
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where ak is the annihilation operator defined in Sec. [TT] 
(extended to d dimensions in an obvious way). 

Introducing an interaction term, H int = 
and the attendant counter term, H c t . = Qai and 

working as in the one-dimensional case, we obtain a shift 
in the poles of the Green function, 


E = 


Ayi 1 

4 W]( 

ker K 


0( A 2 ) 


(C5) 


5 m = —E + 0(A 2 ). For large N, the sum is approximated 
by an integral over the hypercube [0,27r] d . For d = 1, it 
reduces to the previous result, whereas for d > 1 , we 
obtain at lowest order in m and A, 


E = C d \ + ... (C 6 ) 


which is related to the counter-term parameter S m via Numerically, C 2 ~ 0.16, and C 3 « 0.11. 
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